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ABSTRACT 

Recent observations of the Class protostarL1527 IRS have revealed a rotationally sup- 
ported disc with an outer radius of at least 100 au. Measurements of the integrated flux at 870 
^m suggest a disc mass that is too low for gravitational instability to govern angular momen- 
tum transport. However, if parts of the disc are optically thick at sub-mm wavelengths, the 
sub-mm fluxes will underestimate the disc mass, and the disc's actual mass may be substan- 
tially larger, potentially sufficient to be self-gravitating. 

We investigate this possibility using simple self-gravitating disc models. To match the ob- 
served mass accretion rates requires a disc-to-star mass ratio of at least ~ 0.5, which produces 
sub-mm fluxes that are similar to those observed for L1527 IRS in the absence of irradiation 
from the envelope or central star. If in^adiation is significant, then the predicted fluxes exceed 
the observed fluxes by around an order of magnitude. Our model also indicates that the stresses 
produced by the gravitational instability are low enough to prevent disc fragmentation. 

As such, we conclude that observations do not rule out the possibility that the disc around 
LI 527 IRS is self-gravitating, but it is more likely that despite being a very young system, this 
disc may already have left the self-gravitating phase. 

Key words: planets and satellites: formation — Solar system: formation — stars: pre-main- 
sequence — planetary systems — planetary systems: formation — planetary systems: proto- 
planetary discs 



1 INTRODUCTION 

, It is now generally accepted that low-mass stars form 

' through the collapse of co ld, dense molecular cloud cores 

, JTerebev. Shu & CassenI Il984h . These cores, however, typically 

1 have rotation rates that me an that all the mass cannot fall directly 

onto the young protostar (ICaselli et al.l 120021) . A mechanism is 

therefore required to transport this excess angular momentum 

away, allowing mass accretion to take place. 

The standard scenario is that a large fraction of the material 
will fall onto a disc around the young protostar, and that this disc 
will provide the mechanism for transporting the angular momen- 
tum outwards. In many astrophysical discs, this angular momen- 
tum transport could be driven by magnetohydrodynamic (MHD) 
turbulence initiated by the magnetorotational instability (MRI; 
ISalbus & Hawlevlll991r) . At very early times, these cold, dense 
molecular cloud cores are, howeve r, unlikely to be suffici ently ion- 
ized to sustain MHD turbulence JBlaes & Balbuslll994l ). Instead, 
it is thought that these discs could be massive and that disc self- 
gravity could provide the dominant transport mechanism during the 
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earliest stages of star formation (|Toomrelll964l: iLin & Papaloizoul 



ll986l : lLaughlin & BodenheimeJll994l) . 

Although discs around very young (class and class I) pro- 
tostars could be massive, they are difficult to identify because they 
are still heavily embedded in their nascent molecular cloud cores. 
Typically, high-resolution interferometic observations at sub-mm 
and m m wavelengths are used to resolve out the surrounding en- 
velope JRodriguez et al.ll2005l : Eisner et al.ll2005l) . In some cases 
these discs do appear to be quite massive (see, for example, 
lOreaves & Ricelt2011) and refer ences therein). Interferometric ob- 
servations bv lMaurv et al.l ( 1201(1) at 1.3 mm, however, suggest that 
discs around very young protostars are very compact (< 30 AU). 
Magn etohydrodynamic simulations of collapsing molecular cloud 
cores JHermebelle & Tevssien 120081) suggest that even moderate 
magnetic fields can strongly influence angular momentum trans- 
port, and can prevent the disc becoming extended at these early 
times. This may imply that disc self-gravity is not particularly im- 
portant at this stage. 

By contrast, recent observations of the protostar L1527 IRS 
JTobin et al.l 120121) have revealed an edge-on disc with an outer 
radius of at least 100 AU. Since this disc is edge-on, they were 
able to determine the rotational velocity of the material in the disc 
and could confirm that it was indeed in Keplerian rotation about a 
central protostar, with a mass of ~ 0.2 Mq. Not only is this the 
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first confirmation of a disc around a class protostar, it also in- 
dicates that these discs can be reasonably extended (> 100 AU). 
Observations suggest that most of the luminosity is due to accre- 
tion through the disc, giving an accretion rate of 6.6 x 10~^ M© 
yr~^. Continuum observations at 870 /im indicate a disc mass of 
~ 0.007 ± 0.0007 M0, about 30 times smaller than the mass of 
the central protostar. 

If the disc mass estimate is correct, there is insufficient 
mass for this disc to be self-gravitating. Angular momentum 
transport could then be driven by MRI. MHD simulations 
JPapaloizou & Nelsonl200 3) suggest that MRI can pro duce stresses 
resulting in a viscous a CShakura & Sunvaevlll973l) of up to ~ 
X lO"'^. The steady accretion rate in a viscous accretion disc with a 
reasonably low mass and wit h a viscous a of th is magnitude would 
typically be 10"'^ Mq yr"^ JArmitage, Clarke & Palla.2003.) . This 
leaves two possibilities: 

(i) This disc maintains a high accretion rate due to vigorous MRI 
activity supplying angular momentum transport, or 

(ii) The disc maintains a high accretion rate due to the gravi- 
tational instability, and the disc mass is higher than observations 
indicate. 

There is both theoretical and observational evidence 
l lVorobvovil2009l ; iHartmannluOOSb that disc masses around young 
protostars have been systematically underestimated, so option ii) 
may not be particularly surprising. If the disc mass is indeed higher 
than estimated, disc self-gravity could then become a viable mecha- 
nism for transporting angular momentum. As self-gravitating discs 
tend to be m ore centrally condensed than their non -self-gravitating 
counterparts Jciarkel2009HRice & Armitagell2009h . this provides a 
satisfying explanation for their ability to hide mass at high optical 
depths in the hotter inner disc regions. 

In this paper we determine the properties of a self-gravitating 
disc that could provide the observed accretion rate of LI 527 IRS, 
and compare the observational signatures of such a disc with those 
observed. 



2 MODEL 



Here we use t he model developed by Iciarkd ( 120091) (see also 
iRice & Armita ge (2009)). A d isc is susceptib le to the gravitational 
instability if the Q parameter dToomrdl 19641) is close to, but larger 
than, unity. Here we assume that at all radii in the disc 
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(1) 



where Cs is the sound speed, Q, is the angular frequency, G 
is gravitational constant, and E is the disc surface density. 
It is now known that the angular momentum transport in a 
quasi-steady self-gravitating disc can be regarded as viscous-like 
(Balbus &Papaloizou 1999; Lodato & Rice 2004) such that the 
stresses produce an effective viscous a. The mass accretion rate, 
M, in such a quasi-steady disc will satisfy 



M 



SvraCsE 



constant. 



(2) 



The local effective gravitation al a is determi ned by assuming the 
disc is in thermal equilibrium dGammiell200ll) . giving 



(3) 



Using Rosseland mean opacities, k, from ISell & LinI ( 119941) the 
optic al depth can b e estimated as r = kE and the cooling rate is 
then l lHubenvlll990l) 



A = 



ScrT* 
3r 



(4) 



where T is the midplane temperature and a is the Stefan- 
Boltzmann constant. The cooling time is then the thermal energy 
per unit area (CsE/7(7 — 1)) divided by this cooling rate. This is 
then sufficient to determine - given M and stellar mass, M* - the 
sound speed and surface density at any radius in a quasi-steady, 
self-gravitating disc. 

If the disc is subject to irradiation, we modify equation (|4j as 
follows: 



A = 



8cr(T^ - Tit) 
3r 



(5) 



where Tirr represents the temperature of the local radiation field. 



3 RESULTS 

The rotational velocity of the material in the disc around LI 527 IRS 
suggest a central protostar mass of ~ 0.2 M© (Tobin et al. 2012J). 
We use the above model to determine the properties of quasi-steady, 
self-gravitating discs around a central protostar of mass 0.2 Mq. 
We consider two cases: in the first case, there is no significant ex- 
ternal irradiation; in the second, we model the irradiation field from 
the envelope and central star as Tjrr = 30 K, i n keeping with the 
dust temperature estimates of iTobin et al.l ( 120 12h . 



3.1 No Irradiation 

3.1.1 Basic disc properties 

As discussed in Section 2, given M and M*, we can determine 
the surface density and temperature at any radius in a quasi- steady, 
self-gravitating disc. Figure[T]shows contours of disc mass (in Solar 
masses) plotted against accretion rate and outer disc radius. The 
observed accretion rate in L1527 IRS is ~ 6.6 x 10^^ Mq yr^^ 
which, if the outer radius is 100 au, would require a disc with a 
mass of 0.091 Mq. This is shown by the red cross in Figure [T] If 
the outer radius is 150 au, this rises to ~ 0.11 Mq. This is about 
half the mass of the central star and more than 10 times greater than 
that estimated from the observed 870/im flux from L1527 IRS. 



3.1.2 Mass estimates 

We can use our models to determine the 870/im flux for our sim- 
ulated discs. For a system at a distance D, the flux at frequency v 
from an annulus of the disc at radius r from the central star, and 
with radial extent dr is 

2fc 



F^{r)dr 



c^D'^ 



vK{v)T,{r)T{r)2-Krdr, 



(6) 



where k is Boltzmann's constant, c is the speed of light, and T{r) 
is the disc midplane temperature. The frequency dependent opacity 
hi(u) is 



K.{v) — Ho 



(7) 



97(7-l)fcooi«' 



where (3 is typically 1, and Kq is the opacity at a reference wave- 
length of \o — 850/im. 
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Actual Disc Mass 




Figure 1. 2D contours of the actual disc mass calculated from the self- 
gravitating disc models, as a function of the steady state accretion rate M, 
and the disc's outer radius Tout, for a star of mass 0.2A/q. The asterisk 
indicates the required disc mass, assuming an outer radius of 100 au, to 
produce a mass accretion rate like that measured for L1527 IRS. 



Quasi- steady, self-gravitating discs are, how ever, centrally 
condensed JClarkd l200St iRice & Armitagel |2009|) and it is ex- 

S iected that they may b e optically thick even at mm wavelengths 
Greaves & Ricd uOld) . We account for optically thick regions 
(T{r, u) = E(r)K(i/) > 1) by modifying Eq. {6) as follows 



F^{r)dr = 



^v^ K.[v)T.(r)T(r)2Tirdr r sC 1 
^y^I^2TTrdr T > 1. 



i^2 



(8) 



To calculate the integrated flux at 870^m from our simulated discs, 
we simply integrate Eq. ([8} from the inner disc edge (r = 1 au) to 
the outer disc radius, using Ko — 0.035 cm^ g~^ and P — 1. 

The top left panel of Figure|2]shows contours of 870/.im flux in 
mjy plotted agains A/ and outer disc radius, rout. For rout = 100 
au and for M = 6.6 x lO"'^ M© yr"^ the 870/im flux is ~ 475 
mJy. For rout = 150 au, this increases to 650 mJy. The observed 
870/im flux for L1527 IRS is 213. 6± 8.1 mJy. Despite our required 
disc masses (shown in Figure [T} being an order of magnitude - or 
more - greater than that determined for L1527 IRS, the 870/im flux 
is only a factor of 2 to 3 greater than that observed. Our calculation 
is also fairly simple and essentially assumes a face-on geometry. 
The disc in L1527 IRS is edge-on and one might expect a reduction 
in flux at these wavelengths by as much as a fac tor of ~ 2 due to 
this edge-on geometry (e.g., lwhitnev et al.ll2003l) . 

Typically the disc mass is inferred from the long-wavelength 
flux using 



Mdisc = 



D^F^& 



2K,[v)ku'^TA, 



(9) 



Although our models have self-consis t ent di sc temperatures, to 
compare with the results of iTobin et al.l ( l2012h , we use their value 
of Tdust = 30 K. The top right panel of Figure |2] shows the disc 
mass that would be inferred from our simulated 870/im flux. As 
expected, for rout = 100 au and M = 6.6 x 10"'^ M© yr"\ the 
870/im flux would suggest a disc mass of ^ 0.011 M©, almost 10 
times lower than the actual disc mass. Admittedly, if one compares 
the top left and right panels of Figures |2] one will notice that our 
estimate of the disc mass for a flux of 213 mJy is slightly less than 



that determined by 'Tobin et al.] ( 1201 2h (0.00525 M© rather than 
0.007 M©), but this doesn't really change our main result. Using 
the 870/im flux, the mass one would estimate for a self-gravitating 
disc accreting at 6.6 x 10~^ M© yr~^ around a star of M* — 0.2 
M© would be 5 — 10 times lower than the actual disc mass and only 
2 — 3 times higher than that estimated for the disc around LI 527 
IRS.^ 

iTobin et alj ( 120 121) also measure 3.4 mm fluxes which they 
find to be 16.0 ± 1.4 mJy. Assuming /3 = 1 in Eq. Q we can 
calculate the 3.4 mm fluxes from our disc models. This is shown in 
the bottom left panel of Figure |2] For a disc of radius rout = 100 
au and for M = 6.6 x 10~^ M© yr~^ our model predicts a 3.4 
mm flux of ~ 11.1 mJy. This is lower than, but similar to, that 
determined for L1527 IRS. 

We can again estimate the disc mass using Eq. ^ which is 
shown in the bottom right panel of Figure |2] Our simulated 3.4 
mm flux would suggest a disc mass (for properties similar to the of 
L1527 IRS) of 0.018 M©. This is higher than the 870/im flux esti- 
mate (0.011 M©) indicating that, at longer wavelengths, the disc is 
less optically thick. The differen ce is, however, much smaller than 
that seen by iTobin et al.l 12013). Their disc mass estimated from 
the 3.4 mm flux (0.025 ± 0.003 M©) is 3.6 times greater than the 
870/im estimate. They have more confidence in the 870/im flux as 
it is less affected by assumption of /3 but do suggest that this could 
indicate the disc being optically thick at 870/im, but our results 
would suggest that this alone cannot explain the observed differ- 
ence the 870/im and 3.4 mm fluxes. We have assumed 13 = 1, 
reflecting the current uncertainty regarding grain growth, and the 
degeneracy between /3 and the disc mass when interpreting obser- 
vations (Andrews & Williams 2005) . 

We therefore also considered what would happen if we varied 
/3. For P = 0.6 we find that the 3.4 mm flux (at M = 6.6 x lO"'' 
M© yr~^ and rout = 100 au) increases from ~ 11 mJy (when 
/3 = 1) to ~ 16 mJy. If one then assumes /? = 1 when deter- 
mining the disc mass from this flux, this increases the estimated 
disc mass from 0.018 M© to 0.026 M© . The large difference in the 
measured 870/im and 3.4 mm fluxes could, therefore, indicate that 
some grain growth has taken place. If so, this could also suggest 
that such discs would still have optically thick regions even at mm 
wavelengths. 



3.2 Irradiation at Ti,, = 30 K 

In the absence of external irradiation, the gravitational instability 
generates t emperatures in th e outer disc much lower than those es- 
timated by ITobin et al.l ( 12012 ). If we add external irradiation to en- 
sure the outer disc cannot cool below 30 K, how does this affect the 
disc properties and the observable flux and mass? 
Figure |3] shows the actual disc mass as a function of M and 
rout in the presence of irradiation. As irradiation tends to weaken 
the gravitational instability, the disc itself must be more mas- 
sive for a given M — r locus to maintain a marginally unsta- 
ble state dCai et al.ll2008l:lKratter & Murrav-Clavll20lll; iRice et all 



I2OI luForgan & Ricell2013r) . This would push the actual disc mass 
of L1527 IRS to 0.132 Mq, around two-thirds the mass of the cen- 
tral star. The contours are almost vertical, indicating that the disc 
surface density profile is almost independent of accretion rate, and 
depends only on radius. This is suggestive of heating due to the 
gravitational instability being superseded by the heating due to the 
irradiation, and the subsequent weakening of the instability, which 
occurs even at such large disc-to-star mass ratios (IRice & Armitagel 
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Total Disc Flux, A = 870.0 fim 



Estimated Disc Mass, A = 870.0 fim 





Total Disc Flux, A =3.4 mm 



Estimated Disc Mass, A = 3.4 mm 





Figure 2. 2D contours of the self-gravitating disc models observable properties, as a function of the steady-state accretion rate M and the disc's outer radius 
?'out, for a star of mass 0.2Mq, assuming no external irradiation. Top left: the 870 fiva flux (in mjy); top right, the disc mass (in Mq) estimated from the 870 
fiva flux; bottom left, the 3.4 mm flux (in mJy); bottom right, the disc mass (in Mq) estimated from the 3.4 mm flux. Assuming an outer radius of 100 au, the 
red crosses indicate what would be expected for a disc with an accretion rate similai' to that of LI 527 IRS. 



l2009f) . As such, the disc surface density profile is no longer as cen- 
trally condensed as it was in the non-irradiated case. 

Adding this extra mass to the disc's outer regions boosts the 
observed flux by around a factor of 5 (left panels of Figure |4]. As 
the dust mass depends linearly on the flux, the estimated mass is 
also a factor of 5 larger (right panels of Figure 0. 

This is certainly less attuned with the observed flux than the 
non-irradiated case. The uncertainty regarding the temperature pro- 
file of L1527 IRS leaves us uncertain as to the efficacy of self- 
gravitating disc models to explain the observations. While it is 
likely that irradiation does set the temperature of the outer disc, 
modelling of other dust discs in Taurus in the sub-mm suggests 
that the disc temperature at 100 au could be up to a factor of 3 
lower JAndrews & Williams 2005i). 



3.3 Fragmentation 

Although we have primarily considered quasi-steady self- 
gravitating discs, it is thought that if such discs are sufficiently 
unstable they may undergo fragm entation to form bound objects 
dGammiduOOll : iRice et alj|200 3 V This depends on how quickly 
such discs are able to cool and it is expected that it would 
only be possible in the outer parts of such discs (IRafikovluOOSl : 



IStamatellos. Hubber & Whitworthll2007h . The condition for frag- 
ment ation is that the e f fective gravitational vi scosity satisfies a > 
0.06 JRice et alJlJOO^ . iForgan & Ricd ( |201 ih suggest that once a 
exceeds this threshold, fragmentation also requires that the Jeans 
mass changes sufficiently rapidly that a clump would actually form 
on a reasonably short timescale. This allows one to estimate under 
what conditions fragmentation should occur and, if fragmentation 
does occur, the resulting Jeans mass. 



In the non-irradiated case, our calculation show that LI 527 
IRS would need to maintain an accretion rate of around AI ~ 
8 X IO^'^Mq yr~^ or higher to fragment, and is therefore ex- 
tremely close to fragmenting. In the presence of irradiation, this 
critical accretion rate increases, and the disc is therefore stable 
JForgan & Ricel2013h . It is therefore probable that L1527 IRS will 
not fragment but if fragmentation has occuiTed, it may be possi- 
ble to detect enhanced emission consistent with an object of a few 
Jupiter masses or la rger forming via fragmentation in such a disc 
( lGreavesetZll2008h . 



© 0000 RAS, MNRAS 000, 000-000 



A self-gravitating disc in LI 52 7 IRS ? 5 
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Estimated Disc Mass, A = 870.0 fim 
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Figure 4. 2D contours of the self-gravitating disc models observable properties, as a function of the steady-state accretion rate M and the disc's outer radius 
''out, for a star of mass 0.2Mq, assuming irradiation with a uniform temperature of Ti^^ = 30 K. Top left: the 870 /iva flux (in mjy); top right, the disc 
mass (in Mq) estimated from the 870 fiva flux; bottom left, the 3.4 mm flux (in mJy); bottom right, the disc mass (in Mq) estimated from the 3.4 mm flux. 
Assuming an outer radius of 100 au, the red crosses indicate what would be expected for a disc with an accretion rate similar to that of L1527 IRS. 



4 DISCUSSION AND CONCLUSION 

The observation of a rotationally supported disc in the LI 527 IRS 
protostellar system is the first confiimation of a reasonably ex- 
tended (r-Qut > 100 au) disc around a class protostar. The disc 
mass estimated from the sub-mm flux, if correct, suggests that such 
a disc cannot be self-gravitating. 

We show here, that a self-gravitating disc can provide the 
observed mass accretion rate, but the disc mass would need 
to be at least an order of magnitude greater t han that ob- 
served. Such dis c s are very centrally condensed (IClarkel 12003 : 
iRice & Arrrutagd 120091) and so the inner regions are opti- 
cally thick even at long wavelengths. Such discs may produce 
episodic outbursts(Vorobvov & Basu 2005; Vorobvov 2010) and 
the observed FU Orionis outbursts (IHartmann & Kenvotj Il99q) 
may be driven by centrally condensed self-g ravitating discs 
jZhu, Hartmann & Gammi3l2009l ; Izhu et alfcoiof) . 

Using reasonable values for the frequency dependent opac- 
ity, we show that the fluxes produced by the model are similar 
to that observed for LI 527 IRS if external iiTadiation is weak. If 
these fluxes are then used to estimate the disc mass, these would 
give estimates 5 — 10 times lower than the actual disc mass. By 
varying the opacity spectral index /? from its canonical value of 



1 down to 0.6, the fluxes obtained from the model match the ob- 
served values very closely. This might suggest that significant grain 
growth is already underway in the Class phase, perhaps due to the 
self-gravitating disc turbulence reducing the relative g rain veloci- 
ties a l lowin g the grain growth process to be accelerated (IRice et al. I 
I2OO4I ; Iciarke & Lodato ..2009.) . Such a hypothesis remains uncon- 
firmable due to the degeneracies tying j3 to the disc mass, not to 
mention the underlying uncertainties surrounding the grain compo- 
sition and size distribution (Andrews & Williams 2005 and refer- 
ences therein). 

However, in the presence of irradiation, the self-gravitating 
disc model produces fluxes up to 5 times higher than observed. We 
therefore cannot claim that the observations of LI 527 IRS show 
that the disc is massive and self-gravitating, especially as we do 
not model obscuration from the envelope in which the system is 
embedded. If the envelope strongly obscures flux over a wide range 
of wavelengths, then the over-estimated flux produced by the disc 
models in the irradiated case may be greatly reduced, bringing them 
closer to the observed fluxes from L1527 IRS. 

Equally, given the extremely young age of the system and 
the difficulties of modelling the disc in the presence of irradia- 
tion, we can infer that the self- gravitating phase of protostellar discs 
is indeed constrained to the very earliest stages of star formation. 
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Actual Disc Mass 




Figure 3. 2D contours of the actual disc mass calculated from self- 
gravitating disc models in an external irradiation field of Tirr = 30K, as 
a function of the steady state accretion rate M, and the disc's outer radius 
Tout, for a star of mass 0.2A/q. The asterisk indicates the required disc 
mass, assuming an outer radius of 100 au, to produce a mass accretion rate 
like that measured for LI 527 IRS. 



Despite possessing a high accretion rate, the self-gravitating disc 
model's over-estimate of flux in the presence of irradiation sug- 
gests that the disc may already have transitioned from a state in 
which self-gravity dominates as an angular momentum transport 
mode, to a state where MRI now dominates the transport. If this is 
true, and the disc is still somewhat self-obscured, it is evident that 
self-gravity becomes ineffective in protost ellar discs even while th e 
disc-to-star mass ratio remains quite large (iRice & Armitagel2009r) . 
This is illustrated by the disc mass becoming only weakly depen- 
dent on accretion rate in the presence of irradiation, i.e. they are not 
strongly centrally condensed, and are in fact more akin to the typi- 
cal E oc r~^ powerlaw disc profiles We suggest that the actual disc 
mass lies between the lower limit (established by observation) of 
0.007 ± 0.0007Mq and the upper limit (estabhshed by this work) 
of~O.lM0. 
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